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1. INTRODUCTION 

Recently, multilevel inverters (MLI) are widely used in high-power industrial applications because 
of these advantages: low output voltage distortion, reducing the voltage level on semiconductor devices and 
lowering the du/dt changing speed [1]. Especially, in comparison with other MLI topologies such as neutral- 
point clamped [2], flying capacitors [3], multi-modular converters (MMC), the cascaded h-bridge (CHB) [4]- 
[6] has emerged as a prominent one for driving IM system due to its high degree of modularity, allowing 
upgrade and replacement easily. Furthermore, for the applications of induction motor drives, the connection 
of CHB-MLI with multi-phase rectifier transformer ensures the current drawn from the grid have a sinusoidal 
waveform, decreases the total harmonic distortion injected into the grid. 

On the other hand, a considerable drawback of CHB-MLI that produces the common-mode voltage 
(CMV) between the neutral points of the load and the converter [7], which can be calculated by (1) and 
illustrated by Figure 1. Due to this CMV, high-frequency electromagnetic interference (EMI) is introduced, 
causing problems regarding measurement and system monitoring [8]. Furthermore, in IM driving system, the 
CMV generates the common-mode current making detrimental effects on the motor’s productivity and 
endurance [9]. There have been many undergoing researches to minimize the CMV, for instance the addition 
of the EMI active or passive filter [10], [11]. However, this hardware solution can lead to an increase in the 
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size and cost of the system. Therefore, mitigation of CMV based on control method approaches like SVM 
[12], [13], MPC [14], [15] is preferred. 


1 
Vzn = Vemv = 3 (Van + Ven + Ven) (1) 


Among these control methods, multistep MPC has significant benefits since it takes the switching 
nature of the power converter. It is also shown that multistep MPC can improve the steady-state performance 
of the system better than single-step MPC [16]-[20]. However, in [21] multistep model predictive control 
(MPC) with CMV minimization target applies to RL load as the case study only. Regarding IM driving 
system case study, in [22], the multistep MPC controller does not consider CMV minimization target since it 
is designed toward optimizing the switching effort. 

Therefore, this paper proposes a modified multistep MPC compared to [22]. To account for CMV, 
an additional term that penalizes the CMV is added to the cost function before using the SDA algorithm to 
solve the optimization problem. As evidence from the simulation results, this approach also cuts down the 
computational burden, confirms the feasibility of this proposal in practical adoption. 


with xe{A,B,C} 
CHB-MLI ie{1,2...,n} 


Figure 1. IM drive system with CHB-MLI. 


2. MULTISTEP MODEL PRECDICTIVE CONTROL METHOD 

In general, the proposed control method for IM drive system follows the field oriented controller 
(FOC) approach, as depicted in Figure 2 [23]-[25]. There are two major control loops. The classical PI 
controller takes care of the outer speed and flux control loop, whereas the inner current loop is controlled by 
the multistep MPC. Notice that two different weighting factors which account for CMV minimization and 
reduce switching effort are applied to the cost function. 


aS ie Re a Multistep MPC | CHB-MEI 


Switching Lh 
state table 


| Walk) 


Rotor flux model 


Figure 2. The control scheme for induction motor with inner multistep MPC current control loop 
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2.1. Current predictive model 

First, the physical model of the IM drive system is formulated as differential equations, representing 
the relationship between system parameters. After that, a state-space representation of the system can be 
obtained. The state-space representation of the continuous-time prediction model of the IM drive system 
powered by a CHB-MLI can be written as (2): 


faxo = A'(t). x(t) + B'u(t) r 
y(t+1)=C.x(t+1) 2) 
With 
1 1-o 1-0 1-0 
~ (= =e) 0 oTr o w(t) 
1 t= 1-o 1-0 
A(t) = 0 7 (| + =e) w(t) OTr 
1 1 
T. 0 = T, —w(t) 
0 - w(t) = 
2s “1 


f ; i mat. . 

Where the state vector x = [isa isg Yra p,e] € Of, the control input vector u = $ = 
[Sa Sg Sc)’ €U (the set U is the set of voltage level combinations in three phase), output vector y = 
isap = [isa tsp]? € O°, notice that A'(t) is a time-variant matrix because it consists of time-variant 


variables w(t). Specifically, isa ise Vra Vrp stand for stator currents and rotor flux linkages in the 
stationary af -framework. System parameters o, T;,T, and w are defined as total leakage factor, stator time 
constant, rotor time constant and rotor speed respectively. After applying the Forward-Euler discretization 
with a sampling period of T, the discrete-time prediction model can be expressed as (3). 


A=L+T A'E 0% 


n + 1) = A. x(k) + Bu(k) (3) 
B =T.B'e 0% 


y(k + 1) = C.x(k +1) with f 


Based on the prediction model in (4), the control input and output sequence over the prediction horizon N>1 
can be constructed as (5) with relationships in (6). 


Uw) [UK k+ o k+- TU" (4) 

Yw)=[(VK+D)? K+" .. (WFD? | e0™ (5) 
CA CB oo OS 0 

¥(k) = Fx(k) + YUCK) Wherer =|°4"|,y=| Ae e 0 Og 
CAN CA"-1B CAN-2B CAB CB 


2.2. Reference design 

Because the stator current and rotor flux signals are slower than electrical ones, so that over the 
finite prediction horizon N, these parameters are assumed to be constant. The output reference sequence for 
the current controller is denoted by (7). 


VK S [OED OEDS o Ok + NDT € 0? (7) 
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T 
Wherey* (£) =| [i2,(k) + i2,(k). sin( w,€T + po) i2, (k) + i2 (k). sin( w,fT + po + 5 lE 
sp sp 2 


isg(k 
{k + 1;...;k +N}, po = arctan( 28! ) 
isa(k) 


The control input reference sequence (8) is designed while considering a null CMV (Vemy = 0). 


) is the initial angle between i,, and isg. 


U*(k) = [(u* (k) (u*(k+1)" ... (uk +N- 1)" (8) 
2 - 1-0 1—0 i 
. tsJ 2| | disepO 1 1-0) vs Ty o i 
Where u"(Q) = S83 z || = (+E BeO- Lg Prap) 
Cl 3 -v3 sS P —— uw 
= — oO oTr 
2 2 


2.3. Cost function 

To achieve the desired CMV minimization target, the cost function (9) is modified match up to [22]. 
In addition to output reference tracking, the proposed method which is designed to minimize the CMV also 
tracks the input references. 


Iv (ky = Eire 6 + D- yE + DI + AacllAwOll3 + Acmv lu -w (Old) (9) 


In (9),u(€) is the candidate control-input that generates the output current predictiony(€ + 1) = 
igg(€ +1); Au) = u£) — u — 1); acis the weighting factor for lessening switching effort, Acyyis the 
weighting factor for CMV minimization. Representation of (9) in matrix form can be written as (10): 


Iv (k) = IY) — Y°CKD)I + AaclISU(K) — Eu(k — 1113 + cmv llU Ck) — U* CK) (10) 
I, 0, -= 0; i 
-I, l} + 0 a: 
WhereS =| 03 —-I3 + 03 ne te 
ON 05 i 
0; 0; Is disnxan] 3- [3Nx3] 


2.4. Sphere decoding algorithm (SDA) 

In this job, the SDA is adopted to solve the modified optimization problem, to obtain the input 
sequence U(k) = Uopt(k) that minimizes the cost functionJy (k) in detail. The computational process will be 
presented concisely. Further detail about the SDA can be found in [18], [19]. By using (5), the cost function 
(10) becomes: 


Jy(k) = U(k)’ WU(k) +2F(k)'U(k) + €(k) (11) 


Where 
W = Y Y +Ag-STS + Acuvian 
F(k) = Y°Tx(k) — Y Y* (k) — Ag-STEu(k — 1) — AcmyU* (k) 
e(k) = [|| x(k) — Y* (K)||} + AacllEuCk — 1113 + cmv IlU* Ck) I3] 


By solving ol = 0, the optimal solution U„cis yielded in (12), which does not necessary belong to 
the finite control-input set U. 


U„c = —W1F(k) (12) 


Because W is a symmetric and positive definite matrix for Ag,,Acyy>O, there exists only one unique 
invertible lower triangular matrix H, which can be obtained by performing the Cholesky decomposition to 
W-':W-t = HHT. Hence, the matrix H will satisfy the equation: W = HTH. Consequently, the cost 
function (11) can be rewritten as (13). Figure 3 illustrates the graphical representation of a MPC problem in 
two-dimensional space. To simplify, control input vectors referring to the same ellipse will have the same 
cost value. By transforming the cost function (11) into (13), the original ellipses are transformed into circles. 
This procedure plays an important role in SDA. 
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Jv(k) = ||HU(k) — O,,(k)|[Rwhere U,, (k) = HU, (k) (13) 


Figure 3. Representation of the meaning of the matrix H 


The optimal control input for prediction horizon N can be presented as (14): 
Uopt(k) = [Uanc (k) Uane(kK +1) ~ Uare(k +N- 1)" (14) 
Now, reconstruct (14) into a sequence (element by element) as (15): 


Uopt = [U1 u2 + Usw | voi ie {1,2,....3N} (15) 


Figure 4 depicts the evaluation process of each element u;, so-called a node, of the control input 
sequence U in (15). Figure 5 illustrates the flow diagram of SDA. First, an initial sphere Sin; is defined with 
an initial center U,,,(k) calculated by (16) and a selective initial radius Pini. This initial sphere Sini should be 
small enough to contain at least one solution U(k). Then during the evaluation process, the size of the sphere 
is decreased until there is only one solution contained in it, and that is the optimal control input sequence. 
Each time a node is visited, the radius p; is calculated by (17) before evaluating the condition (18). If p; 
violates condition (18), the node u; will be discarded and its following nodes (from u;,,to u3y) will be 
discarded without performing any computation to avoid unnecessary computation. The control input having 
the smallest radius will be the optimal solution U(k) = Uopt(k). 


Popt = Pini 


u] 
u2 
U3N-1 eee eee eee eee an 
U3N eee eee eee eco -N n 


Update Usp »Popt 


(@: Node is visited, &: Node is discarded, ®: Node is discarded without performing any computation) 


Figure 4. A tree-diagram of the control input sequence 
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Uaclk), H, Ppast =Pinis 
i=0, co=0, po=0, u(k-1) 


2 i= | H hf 
pi= (Hj rU priyUuctij) +Pi-1 


i=i-1,co=cotl 


N 
apply Uodk) 


Figure 5. SDA flow diagram [20] 


Ouc(k) = H(-W*)F(k) (16) 
pr = || (hia... hiu] [uy... uj)” — ŭuci ||3 + pfa (17) 
ee N ` I 
Hļ{i1:i] Uopt(1:i] Ductil 
Condition: p? < Popt (18) 


2.5. Switching state table 


When the optimal solution U(k)is obtained, only the first element u,,,-(k) is employed to the find 
the switching states for the inverter. Notice that u,,,(K) is the phase voltage levels. Table 1 shows the 
switching states selection strategy clearly. In Table 1, variable represent the phases, the index notation 
denotes the position-number of a H-bridge in phase, is the voltage level of phase, is the voltage level of the 
H-bridge in phase, are the switching states of valve 1 and 3 of the H-bridge in phase. Figure 1 illustrates this 
switching state selection strategy. 


Table 1. Switching state selection strategy 
S Sxi (Sas Saia) 
X Sx Sag TN Sy2 (Soi Sa) Sx3 Bax Saa) Sx4 (Sau Saa) Sys (Seii Sea) 
:0) 1; ; 0) 1; 


+5 
+4 1(1;0) 1(1;0) 1(1;0) 1(1;0) 0(0;0) 
+3 1(1;0) 1(1;0) 1(1;0) 0(0;0) 0(0;0) 
+2 1(1;0) 1(1;0) 0(0;0) 0(0;0) 0(0;0) 
+1 1(1;0) 0(0;0) 0(0;0) 0(0;0) 0(0;0) 
0 0(0;0) 0(0;0) 0(0;0) 0(0;0) 0(0;0) 
-1 0(0;0) 0(0;0) 0(0;0) 0(0;0) -1(0;1) 
5) 0(0;0) 0(0;0) 0(0;0) -1(0;1) -1(0;1) 
-3 0(0;0) 0(0;0) -1(0;1) -1(0;1) -1(0;1) 
-4 0(0;0) -1(0;1) -1(0;1) -1(0;1) -1(0;1) 
5 -1(0;1) -1(0;1) -1(0;1) -1(0;1) -1(0;1) 
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3. SIMULATION RESULTS 

The proposed control scheme is simulated by MATLAB/Simulink for IM drive system with an 11- 
level CHB, DC voltage supply per CHB Vg, = 600V and the prediction horizon N=3. The IM drive system 
parameters and the simulation scenarios are shown in Tables 2 and 3 respectively. To definitely demonstrate 
the benefits of the proposed control method in company with [22], the CMV weighting factors are applied to 
the cost function from 1,35s to 1,45s. The sampling period T = 50us. 


Table 2. IM drive system parameters 


Nominal Power Pam 1119 KW 
Nominal Electromagnetic Torque Mam 6906 Nm 
Nominal Current lam 234.18 A 
Nominal Voltage Uam 3.3 kV 
Nominal Speed Wam 1470 rpm 


Table 3. Simulation scenarios 


Time 0-0.5 [s] 0.5-1 [s] 1-1.25 [s] 1.25-1.35 [s] 1.35-1.45 [s] 
Reference speed Magnetization Gradually increase speed Speed Stabilization W,.¢ = 1470rpm 
Wrop = 0 Wref = 1470rpm 
Weighting Factors Aac ~ 0; Acuy ~ 0 Àcmy ~ Oac =4_ Acmy = 4Àac = 4 


It is presented in Figure 6 that the rotor flux, the rotor speed, and the electromagnetic torque are 
smoothly ramped-up and quickly reach their reference values. Figure 7 illustrates the stator current and the 
CMV at the steady-state. The effectiveness of the proposed control method is clearly demonstrated when the 
CMV weighting factor is applied in the period of (1.35-1.45) [s]. Especially in (1.15-1.25) [s] and (1.25-1.35) 
[s], the CMV is registered at 1800V (3V4c), while in (1.35-1.45) [s], the CMV is drastically reduced to 200V 
G. Vac). Moreover, Table 4 depicts the number of evaluated nodes during the SDA optimization process. 


It can be seen that this number is significantly decreased from 144,969 nodes to 88,396 nodes, 
indicates the reduction of the computational burden evidently. However, these benefits are traded for a small 
increase of THD and switching effort, specifically from 0.58% to 0.61% for THD and from 2.43 to 2.78 
times per cycle for switching effort. Finally, Figure 8 shows that the stator current in ap -framework is 


attained its reference for all simulation scenarios with the error below 5%. Hence, the additional CMV 
minimization target does not affect the system performance. 


0.2 0.4 0.6 0.8 
Time (s) 
GBB Reference value J Feedback value 


Figure 6. Graphs of motor speed, torque and motor flux 


Modified multistep model predictive control for three-phase induction motor drive system ... (Bao Binh Pho) 


2258 O ISSN: 2088-8694 


MADAAADAANARAARRAALRAALDDALRAAI AAI LAL AL LALLA LIALNIVINVAS: 
AAA DATAAA AAA RIAA AAA A AAA A AN AA AR AAAAARAAAAA RAI: 
IVVVVVVVVVVVVVVVVVVYVVVVVVVVVVVVVV VV VV VV VV V VV VY: 
LN NN NN NN : 


CMV (V) Current (A) 
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Time(s) 
E Phase A E PhaseB = Phasec M Common-mode Voltage 


Figure 7. Graphs of stator currents and common-mode voltage 


Current aB (A) 


Error (%) 


Time (s) 
WB iso ref Wisg ref Iisa feedback (is feedback [Error Isa QJ Error isp 


Figure 8. Graph of stator current and error value in the aß- framework 


Table 4. Results about THD, switching effort, and number of evaluated nodes 


Time 1.15-1.25 [s] 1.25-1.35 [s] (Reference [22]) 1.35-1.45 [s] (Proposal) 
THD 0.46% 0.58% 0.61% 
Switching Effort 3.26 (times/ cycle) 2.43 (times/ cycle) 2.78 (times/ cycle) 
Number of evaluated nodes 314336 (nodes) 144969 (nodes) 88396 (nodes) 


4. CONCLUSION 

This paper proposes a modified multistep MPC for IM drive system with a CHB-MLI. This 
proposed method achieves two main contributions: firstly, the reduction of CMV while retaining the current 
quality and switching effort at the desired range. Secondly, the required execution time for microcontrollers 
is decreased. Based on this potentiality, future research work will be focused on improving the controller 
performance in advance for practical implementation. 
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